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Abstract 



Q ' The temperature dependence of the rate of the reaction CH4+H -^ CH3+H2 is 

c/3 ■ studied using classical collision theory with a temperature-dependent effective 

• 1— I . potential derived from a path integral analysis. Analytical expressions are ob- 

J>^i tained for the effective potential and for the rate constant. The rate constant 

i-^ \ expressions use a temperature-dependent activation energy. They give better 

agreement with the available experimental results than do previous empirical 

fits. Since all but one of the parameters in the present expressions are ob- 

^ . tained from theory, rather than by fitting to experimental reaction rates, the 

expressions can be expected to be more dependable than purely empirical ex- 

(^ \ pressions at temperatures above 2000 K or below 350 K, where experimental 

^H ' results are not available. 
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I. INTRODUCTION 

The validity of the Arrhenius expression for the rate constant of a bimolecular reaction, 
k = Aexp{—Ea/RT), where Ea is the activation energy and R is the gas constant, has long 
been the subject of controversy ||1[. -Ea is generally assumed to be independent of temper- 
ature, and any deviation from linear behavior in the plot of log A; vs. 1/T is attributed 
to temperature dependence in A. Classical collision theory 0] suggests only a weak tem- 
perature dependence of the preexponential factor, of the form A oc T^'^, while theoretical 
arguments based on transition-state theory (TST) with corrections for quantum mechanical 
tunneling |l|J^ can predict a significantly different dependence. 

We consider here the reaction 

CH4 + H ^ CH3 + H2. (1) 

There exists an extensive literature of theoretical and experimental studies 0-14] for this 
reaction on account of its importance as an elementary reaction in hydrocarbon pyrolysis. 
The experimental results span the temperature range from 372 K up to almost 2000 K. 
TST studies [||-0 have predicted a distinct upward curvature in the Arrhenius plot, and 
an analysis by Shaw p[ of experimental results through the year 1978 generally supported 



this prediction. Subsequently, Sepehrad et al. |12] concluded that in fact a linear Arrhenius 



plot was more consistent with the data, after omitting some apparently unreliable earlier 
results and adding new results of their own. More recent analyses by Baulch et al. [0 and 



by Rabinowitz et al. [|T^ discerned curvature with A oc T^ and A oc T^'^^, respectively. 

We propose an alternative theoretical approach to this problem, based not on transition- 
state theory but on classical collision theory with a temperature-dependent effective potential 
energy function. This approach is based on a path-integral analysis developed by Feynman 



15| , p^ for problems in quantum statistical mechanics. He noted that the equation for the 
statistical density matrix is formally identical to that for the kernel that expresses the time 
dependence of the wavefunction of a quantum mechanical particle over a time interval that 
is taken to be negative and imaginary. Thus, calculations in statistical mechanics can be 
carried out using the path integral techniques of quantum dynamics. This idea can be used 
to derive an effective potential \4fj for the chemical reaction such that a classical mechanical 
calculation of the reaction rate with the effective potential is approximately equivalent to 
a quantum mechanical treatment using the actual potential V fl^. This simple approach 
was used previously to describe the diffusion of H on a Cu surface |]TB[ , and the results were 
later found to agree with those from an elaborate reaction-path variational TST calculation 
with semiclassical adiabatic tunneling corrections |I^]. 

The path integral analysis transforms V into a temperature- dependent iunction Ves{T), 
from which we obtain a temperature-dependent activation energy. Thus, our model for the 
chemical reaction is classical collision theory but with a temperature-dependent Ea. We 
will include no temperature dependence in A other than the classical T^/^ factor. We will 
introduce a slight modification into the analysis so that the path- integral result, which is 
derived as a perturbation theory about the high-T limit, can be smoothly interpolated to 
the correct low-T result. For V we will use the empirical hydrocarbon potential of Brenner 

20| , which is constructed from Morse-type functions with modifications to take into account 



non-local effects. 



II. METHOD 

To calculate the partition function from the statistical density matrix it is sufficient to 
consider only paths that return to their starting point. The path integrals that need to be 
evaluated are very difficult on account of the large number of degrees of freedom needed to 
describe the many possible paths. However, if h/kBT is small, where ks is the Boltzmann 
constant, then one can derive a simple approximate expression for the partition function, 
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-PV.s{'-c)dx, (2) 



where P = l/fc^T and 



' inksT' 



\/2na J-Qo 



a^ = nyUmkBT. (4) 

The significance of this result is the fact that Eq. (|^) has exactly the form of the classical 
partition function except that V is replaced by the effective potential V^s, which is just a 
Gaussian average of V with a temperature-dependent standard deviation a. Equation (|[) 
is for a system with only one degree of freedom, but the extension to an arbitrary number 
of degrees of freedom is straightforward. 

It has been suggested []T^ that classical dynamics on the potential Ves could be used 



to simulate quantum dynamics on the true potential V. This would not be valid for the 
dynamics of a single particle, but it is a reasonable hypothesis for a statistical ensemble of 
particles as in a molecular dynamics simulation of a chemical reaction. In fact, it can be 
shown [|l^ that when Ves is used in place of V in TST one obtains the standard Wigner 
tunneling correction. 

This approach provides an appealing qualitative model for quantum effects. At a mini- 
mum of V, averaging over neighboring points according to Eq. (j^) will increase the potential. 
This accounts for the fact that the minimum of V is in practice inaccessible to the system 
on account of the impossibility of localizing a quantum mechanical particle. In effect, the 
averaging provides a zero-point energy correction. At a maximum, the averaging reduces V. 
In effect, this is a tunneling correction. At a saddle point of V, averaging over a given coor- 
dinate will reduce the potential if V" is negative and increase it if V" is positive. Elsewhere, 
V can be approximated as a linear function, in which case the averaging in Eq. (§) will 
have little effect. At high T the quantum effects will be small, because if the average kinetic 
energy is large then it is only rarely that a particle will be close enough to the potential 
surface to sense the difference between V and Kg. Accordingly, the standard deviation given 
by Eq. (^) goes to zero in the limit of infinite T. 

At low T the approximations used to derive Eq. (^ can be expected to lead to a significant 
error, since the system will spend much of its time in the quantum regions of V. However, a 
minor modification of the theory will ensure that it give the correct low-temperature limit. 



At T = K the system will be at rest at the nearest local minimuni of V. The energy of 
the system will be the value of V plus a zero-point energy correction. A minimum of V in 
principle corresponds to a stable chemical species, for which the zero-point energy can be 
determined empirically from analysis of the vibrational spectrum. Thus, we can replace T 
in Eq. (§) with 
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T + To, (5) 



where Tq is a constant chosen such that Eq. (j^) reproduces the empirical zero-point energy 
of a known species. 

Now consider the application of this theory to the hydrogen abstraction reaction, Eq. (p. 
We will treat this as a problem in two degrees of freedom {x,y), where x is the C — H 
distance for the reacting hydrogen and y is the H — H distance for the reacting hydrogen. 
We will assume that the C — H — H configuration is linear along the reaction path. The 
angle 6 between the reactive and nonreactive C — H bonds will be treated as a quadratic 
polynomial in x that interpolates between 109.5° for CH4, 120° for CH3, and the transition- 
state geometry of 102.4° at x = 1.08 A [0. Then, 

Ves{x,y) = -^— fdzye~'y/'''y f dz^e-'"/^''"V{x + z^,y + Zy), (6) 
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where 
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Note that the mass in Eq. (H) has been replaced by the reduced mass fi^ = mn^l+mn/mc)'^ 
or fiy = niH/^. 



For V{x, y) we use Brenner's potential I pO|) which is a sum of two-body interactions 



V{x, y) = Vcr{x, y) + VnnijJ, x). (8) 

The VcH and Vhh each have the form 

V^r, q) = /.(r)|Dr^ exp [ - (3^^ (r - rf^)] - B,{r, q)D\''^ exp [ - (3^^^ (r - rf^)] }, (9) 

where r is the coordinate of the primary, two-body, interaction while q is the coordinate of 
the "environment." Equation (^ has the general form of a Morse potential but is modified 
by the functions /, and Bi. fi{r) is a cutoff function that smoothly interpolates to zero in the 
limit of large r. Bi{r,q) is a rather complicated function that models the effects of nearby 
atoms on the primary interaction. For given (x, y) Eq. (|^) can be evaluated by numerical 
quadrature. The values 

To,, = 92i^, To,, = 582K (10) 

for the temperature shift parameters give agreement with the spectroscopically determined 
zero-point energies for CH4 [^ and H2 [^, respectively. 



There is a minor inconsistency in using this procedure with the Brenner potential. For 
the C — H well depth Brenner simply used the bond energy of the CH molecule without 
subtracting the zero-point energy. Thus, we are in a sense adding zero-point energy to a 
potential that already includes it. This will make the calculated activation energy smaller 
than it ought to be. We will assume that any errors from this procedure will be insignificant, 
since the quantum effects for C — H interactions are much smaller than those for H — B. 
interactions. However, in principle, one ought to refit the potential using the correct well 
depth. For the H — H potential, which causes most of the quantum effects, Brenner did use 
the correct well depth for H2. 

We find that it is possible to accurately approximate Ves with an analytic expression. 
Let Vcs = Vcn + Van, where 
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Note that Eq. (|1^) can be written as 
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where 
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The term V^ simply ignores the coordinate dependence of fi and Bi for purposes of evalu- 
ating the integral. The integrals in Eq. (|TB|) can be evaluated exactly, giving an expression 
identical to that for Vi in Eq. (|]) except with the prefactors Di replaced by the temperature- 
dependent parameters 
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For the C — H interaction, we will use Vch ~ V^h ^^^ ignore the three correction terms. 
For the H — H interaction, since a^ is over twice as large as cr^, we will include all of the 
correction terms but replace them with approximate analytical expressions. Let 



r(l) 
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V^\i = <j,'j^>{y) + BMy,x)<l>'^>{y) 



W, 



(18) 



where 
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4^)(y) = Di^ijdz [Mniy + ^) - hniy)] exp [-(3^''\y 

i«(y) = -DJjU / dz [funiy + z)- fnuiy)] exp [-/?(^) (y - rj^ + z) - zV2a, 



(19) 
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We find in practice that the functions 0)j and (f)\ can be fit quite accurately as sums of 
Gaussians, 



^w ^ /,{i.i)e-(i/-cL^'^0'M^'^'' + ;^a,2)g-(.-c(^'^))'M^-^)'^ 



(21) 



with the parameters h, c and if given by constants or by polynomials in (T + To^y)~^. The 
coefficients of these polynomials can be obtained by fitting to a set of exact values of the 
(p^^^ from numerical quadrature of Eqs. (|T9|) and (pOl). 

Another function that can be fit in terms of a sum of Gaussians is 



0(2)(y,x) = -4Ue-^'"'(^-4^^)|d^ [i?HH(2/,x + z) -5HH(y,x)] e-^'/2- 
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'^(2,l)g-(y-c(2,i))'/«;(2,i)^ ^ ^(2,2)g-(3;-c(2.2))V^«(2.2) = 



(22a) 
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In this case the parameters are polynomials in (T + T^y) ^ and y. The second correction 
term is given by 
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Finally, we express the third correction term as 

V^'\y.x) = fcuix)fnniy)ci^^'\x,y)/iT + To,y), 



(23) 



(24) 



with 0^^) fit as a polynomial in x and y. 

Using polynomials of at most degree 2 for the parameters we can fit the exact Gaussian 
averaged potential with an accuracy of ±0.005 kcal/mol in the vicinity of the transition 
state for T as low as 300 K. A table of fitting parameters and a full error analysis for the 
the (f)^^^ will be presented elsewhere. 

III. RESULTS 

We obtain a temperature-dependent activation energy as the difference between Ves 
evaluated at the CH5 saddle point and at the reactants well. In Fig. 1 the sohd curve shows 
our quantum mechanical result for Ea, calculated from the Gaussian average with standard 
deviations ai in terms of effective temperatures according to Eqs. (|^. This result is in effect 
an interpolation between the high-temperature path integral analysis and the empirical low- 
temperature limit. For T > 300 K we find that Ea can be accurately fit with a quadratic 
polynomial in T~^, 
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Ea{T) ^ 12.07 kcal mol"^ - (741.2 kcal mol"^ K) T"^ + (7.47 x 10^ kcal mor^ K^) T'^. 

(25) 

The dashed curve shows the result of a purely high-temperature analysis, with the Uj in terms 
of the actual temperature. The two results are in general agreement down to approximately 
600 K. (The unsteadiness in the curves at low T is due to uncertainty in the determination 
of the geometry of the activated complex.) 

The corresponding Arrhenius curves are shown in Fig. 2, where they are compared with 
classical collision theory (the dotted curve, corresponding to a temperature-independent 
activation energy equal to liinT-,00 Ea{T) ), with a recent multidimensional semi classical 
variational TST analysis 0] (the dash-dot curve), and with experimentally determined 
reaction rates [^|T2|,|T1|. Our theoretical treatment yields the activation energy, but not the 
Arrhenius prefactor, which must be fit to the experimental results. The theory for the T 
dependence of Ea is exact in the high-temperature limit. Therefore, it is best to consider 
only high-temperature rate constants for the purpose of determining the prefactor. We will 



assume that the measurements by Rabinowitz et al. ||Tj], from the recent Brookhaven flash 
photolysis-shock tube study, are the most reliable of the high temperature results. Fitting 
the classical collision theory expression 

k = AoT^/^e-^"/^^ (26) 

to the the Brookhaven results, with our theoretical formula for the function Ea{T), gives 
log]^Q(Ao /cm^s~^) = —11.536. It is clear from Fig. 2 that the Gaussian average with Tes is 
in better agreement with the low temperature measurements than are the Gaussian average 
with the actual T or the classical theory with Ea independent of temperature. 



IV. DISCUSSION 

A. Comparison with transition-state theory 

It can be seen in Fig. 2 that our theory (the Gaussian average with Tgfj) and Truong's 
TST [|10| both predict a distinct upward curvature in the Arrhenius plot at low temperature. 



They disagree at high temperature, with our analysis predicting a slower reaction rate than 
that predicted by the TST study. 

Our theoretical approach differs from TST in the way in which we describe the classical 
dynamics and in the way we include quantum effects. In TST the expression for the rate 
constant can be expressed approximately in the form k = aoT^Te~^'^^^^ , where Oq and 
X are temperature-independent constants and T{T) is a quantum correction |§|,^. The 
reactants and the transition state are assumed to be in equilibrium, and the factor T^ comes 
from the resulting equilibrium constant, expressed as a ratio of partition functions. For 
CH4+II ^ CH5 the value of x is 3/2. Thus, we can express the TST result in the form 

k = aoT^/^e-^'"/^^, (27) 

with Ea equal to a constant for classical TST and Ea{T) = -E^^'^ In (r(T)) for quantum TST. 
Our approach is to use collision theory for the classical dynamics, which gives x = 1/2, 



and to determine Ea(T) directly from the temperature-dependent effective potential instead 
using one of the standard semiclassical expressions for r(T). However, we could just as well 
have used the effective-potential values for Ea in Eq. ( pTj) instead of using the collision-theory 
expression Eq. (^61) . 

A clearer comparison of theory and experiment is given by Fig. 3, which shows the 
difference between logj^g ^ ^^^ log^g ^ci, where kd is the rate constant predicted by classical 
theory. We obtain kd from Eq. ( [2B| ) by setting Ea = 12.07 kcal/mol, which is the value 
given by the Brenner potential. We use the same value of Aq as with the T-dependent 
Ea, since the quantum and classical theories have the same high-T limit. The solid curve 
in Fig. 3 corresponds to our quantum collision theory, using Eq. p6| ) with the effective- 
potential Ea{T) and with Tcs in the standard deviations for the Gaussian averaging. The 
dotted curve results from using this same Ea{T) in the TST expression Eq. (pTj). The solid 
curve is clearly in better agreement with the experimental points. 

Figure 3 also shows Truong's variational TST results [T^ (the dash-dot curve), which 



were based on ab initio calculations of the potential, and it shows results from a varia- 
tional TST study by Joseph et al. ^ (the dash-dot-dot curve), which used an analytical 
semiempirical potential energy function. The prefactor oq used for our effective-potential 
TST was chosen to give agreement with the high-temperature experimental results. Note, 
however, that different values of Oq would lead to qualitative agreement between the effective- 
potential TST and the variational TST studies. Although there is a fair amount of scatter 
in the experimental results, our quantum collision theory and the TST studies appear to 
agree equally well with the low-temperature experimental points, where quantum effects are 
expected to be most important. This is consistent with previous studies for other systems 
|l^,19|, which showed agreement between quantum effects calculated from semiclassical TST 



quantum corrections and those calculated from effective potentials. 

At high temperatures the results from Truong's TST results lie above most of the experi- 
mental points. The results from semiempirical TST are in better agreement with experiment, 
but extrapolation of them to higher T would give a larger rate constant than does our col- 
lision theory. Since barrier recrossing becomes more likely at higher T, one can expect that 
TST will increasingly overestimate the reaction rate as T increases. Truong used an ah initio 
result for the barrier height (16.3 kcal/mol) that was higher than the value used by Joseph 
et al. (12.9 kcal/mol) or the value used by us (12.07 kcal/mol). However, as T increases, 
the value used for Ea becomes less significant than the value used for the exponent x. At 
high T the assumption on which collision theory is based, that only the translational degrees 
of freedom of the colliding molecules need to be taken into account, is reasonable since the 
lifetime of the activated complex will typically be too short for there to be significant con- 
version of translational energy into vibrational or rotational energy before it breaks apart. 
At lower temperature TST will be more accurate than collision theory, but the error intro- 
duced by using a; = 1/2 will become relatively less important than errors in the value of Ea- 
The excellent agreement with experiment from our quantum collision theory would seem to 
indicate that this is what happens for the particular reaction considered here. In any case, 
the approach we have presented for calculating quantum effects can be used with TST as 
well as with collision theory, and it is conceivable that for other reactions, at intermediate 
temperatures, TST might be the more appropriate choice. 



B. Comparison of analytical expressions for the rate constant 

In Fig. 4 we compare our calculated results with various "best fits" of the experimental 
points P.|T^-[14|. The fitting functions are given in Table 1. For Fig. 4 we computed our 



rate constants using analytic expresions for Ea(T). The solid curve was computed using 
the quadratic fitting function in Eq. (pS)). This curve is almost indistinguishable from the 
corresponding curve in Fig. 3, which was computed from our actual values for the Ea- The 
dashed curve results from a linear fit for Ea{T), which yields an expression for k that has 
only 3 parameters. Most of the empirical fits in Fig. 4 are also 3-parameter expressions. 

Our results in Fig. 4 are in better agreement with the experimental points than are any 
of the empirical expressions. The expressions of Shaw P] and of Baulch et al. |jT^ (the dash- 
dot and dash-dot-dot curves, respectively) seem to be too low in the low-T region. Both 
of these fits were carried out before the data from the high-T measurements by Rabinowitz 
et al. \Ti\ were available. Shaw's fit agrees with these new data while the fit by Baulch et 



al. appears to be too high in the high-T region. Not surprisingly, the fit by Rabinowitz 
et al. (the dash-dot-dot-dot curve) seems accurate at high-T. However, it appears to be 
too high at low T. Note that only one of our parameters (the prefactor Aq) is fit to the 
experimental results. The remaining parameters are calculated theoretically. (Of course, 
other experimental results, such as the bond dissociation energies of CH4 and H2, enter our 
theory indirectly since we based our calculations on an empirical potential function ||20|| .) 

The most significant difference between our expressions for k and the various fits is 
the behavior that is predicted at very high T, above 2000 K. Our expressions predict a 
significantly lower reaction rate in this region than do any of the others. Unfortunately, 
there are no available experimental results above 2000 K with which to compare. The 
quantum effects at very high T are insignificant. Therefore, the accuracy of our results in 
that region depends only on the accuracy of the activation energy given by the Brenner 
potential and on the accuracy of the Arrhenius expression from classical collision theory. 

C. Conclusions 

We have shown that the Arrhenius expression from classical collision theory with Ea 
expressed as a polynomial in 1/T gives agreement with the experimental measurements 
for the rate constant for CH4+H that is at least as good as that from fits that include T 
dependence in the prefactor. This is especially striking since only one of the parameters in 
our expressions is determined by fitting to the experimental points. The fact that the other 
parameters are determined theoretically suggests that our expressions are more likely to be 
dependable at temperatures above 2000 K or below 350 K, where experimental results are 
unavailable. 

These calculations require that the potential energy surface be known. Otherwise, there 
is no choice but to choose a functional form for k{T) containing parameters that are all 
determined by fitting to measured reaction rates. Traditionally, an Arrhenius expression is 
used with the temperature dependence in the prefactor. However, the theory developed here 
provides a justification for choosing a functional form with temperature dependence in the 
activation energy. 



The success of our theory at describing this reaction rate over a wide temperature range 
has two important imphcations for molecular dynamics simulations of chemical reactions. 
First, it supports the use I^S] of the Brenner hydrocarbon potential in those studies. This is 
an empirical potential, parameterized to a data set consisting of properties of stable chemical 
species, yet it seems to be able to accurately model the unstable CH5 transition state, since 
the curvature of the Arrhenius plot depends sensitively on the topography of the potential 
in the vicinity of the saddle point. 

The second implication is that it is possible to use a classical molecular dynamics compu- 
tation, with a Gaussian-averaged effective potential, to model processes in which quantum 
effects are important. The evaluation of the potential function is the most time-consuming 
step in these computations. Therefore, the use of numerical quadrature to perform the 
Gaussian averaging would be impractical. However, the use of our analytic approximation, 
given by Eqs. ([I2|) through (|2^) , will have only a minor effect on the computational cost. 
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TABLES 



TABLE I. Expressions for the rate constant for CH4+H— > CH3+H2. 



rate constant, cm'^ molecule s 



1 0-1 



method 



reference 



k = 3.02 X 10-1° exp(-6631/T) 

k = 2.35 X 10-i'^r2exp(-4449/r) 
k = 2.18 X 10-20 r3exp(-4045/r) 

/c = 6.4 X 10-18 7^2.11 exp(-3900/T) 

k = 2.63 X 10-12 ri/2 exp ( - 6076/r 

+3.730 X IOVT^ - 3.76 x lO'^/T^) 



empirical fit 



empirical fit 



Sepehrad 
et al."^ 



Shaw^ 



empirical fit 




Baulch 

et al." 


empirical fit 




Rabinowitz 
et al^ 


cohision theory with 




this study^ 


effective potential; 






quadratic expression 


ioi Ea 





k = 2.63 X 10-12 ri/2 exp ( - 6076/r 

+2.759 X 10Vt2) 



collision theory with 
effective potential; 
linear expression for Ea 



this study'^ 



-Ref. @. 
^Ref. il. 



^Ref. [|13|. 
•iRef. mi. 
''Prom Gaussian average of potential with T^s in the standard deviations. 
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FIGURES 

FIG. 1. Temperature dependence of the activation energy. The dotted line shows the classical, 
temperature-independent, result. The dashed curve was calculated from the Gaussian-averaged 
potential using the standard deviations given by Eq. (Q) . The solid curve was calculated from the 
Gaussian-averaged potential with T in Eq. (^ replaced by T^s according to Eqs. ^ and (p^). 



FIG. 2. Arrhenius plot of rate constant vs. 1000 K/T for the reaction CH4-|-H^ CH3-I-H2, 
with k in units of cm^ molecule"^ s~^. Experimental measurements, from Refs. |5|, |12], and |14|, 
are indicated by dots. The dotted line corresponds to classical collision theory, while the dashed 
and solid curves correspond to our quantum theory using T and Tcfj, respectively, in the standard 
deviations. The dash-dot curve shows results from Truong's variational TST In 



FIG. 3. Difference between log^o k and log^Q kf.\ vs. 1000 K/T, where fed is the rate con- 
stant from classical collision theory with temperature-independent Ea- The units of k are 
cm^ molecule"^ s~^. The symbols indicate experimental points from the following references: Ra- 
binowitz et al. [Q (-I-); Sepehrad et al. |l^ (□); and the various studies reviewed by Shaw |^ 
(o). The dashed and solid curve shows results from the present study, using Tgfj in the standard 
deviations of the Gaussian average. The dotted curve corresponds to classical TST, according to 
Eq. (27), with Ea{T) from the Gaussian average. The dash-dot curve corresponds to Truong's 



variational TST |jlO|, which uses ab initio calculations for the potential energy. The dash-dot-dot 
curve corresponds to the variational TST of Joseph et al. ^, which uses a semiempirical potential 
function ("J3"). 

FIG. 4. Difference between log^Q k and log^Q k^i vs. 1000 K/T, using the analytical expres- 
sions for k{T) in Table 1. k^i is the rate constant from classical collision theory with tempera- 
ture-independent Ea- The units of k are cm^ molecule"^ s~^. The symbols indicate experimental 
points, labeled as in Fig. 3. The dotted, dash-dot, dash-dot-dot, and dash-dot-dot-dot curves are 
the empirical fits from Sepehrad et al. [12|, Shaw [^, Baulch et al. ||l^, and Rabinowitz et al. |14|, 



respectively. The solid and dashed curves are from the present study with a quadratic and a linear 
fit, respectively, for Ea. 



13 



